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ABSTRACT 

We argue that combined observations of galaxy rotation curves and gravitational lens¬ 
ing not only allow the deduction of a galaxy’s mass profile, but also yield information 
about the pressure in the galactic fluid. We quantify this statement by enhancing 
the standard formalism for rotation curve and lensing measurements to a first post- 
Newtonian approximation. This enhanced formalism is compatible with currently em¬ 
ployed and established data analysis techniques, and can in principle be used to rein¬ 
terpret existing data in a more general context. The resulting density and pressure 
profiles from this new approach can be used to constrain the equation of state of the 
galactic fluid, and therefore might shed new light on the persistent question of the 
nature of dark matter. 

Key words: equation of state - gravitational lensing - methods: data analysis - 
galaxies: halos - galaxies: kinematics and dynamics - dark matter. 


1 INTRODUCTION 

One of the most compelling issues of modern astrophysics is 
the open question concerning the nature of the dark matter 
which dominates the gravitational field of individual galax¬ 
ies and galaxy clusters. [See for instance Persic et al. (1996), 
Borriello & Salucci (2001), and Salucci & Borriello (2003).] 
While the current consensus in the astrophysics community 
is to advocate the cold dark matter (CDM) paradigm, no di¬ 
rect observations of the equation of state have been carried 
out to confirm this widely adopted assumption. Efforts to 
confirm this assumption include attempts to detect elemen¬ 
tary particles that have been suggested as cold dark matter 
candidates. However, experiments that aim (for instance) to 
detect massive axions with Earth-based detectors (Particle 
Data Group: Eidelman et al. 2004, §22.2.2) do not yet yield 
a positive result. 

A different approach to analysing the nature of dark 
matter has been suggested by Bharadwaj & Kar (2003) 
who first proposed that combined measurements of rotation 
curves and gravitational lensing could be used to determine 
the equation of state of the galactic fluid. Whereas their 
analysis made particular assumptions on the form of the ro¬ 
tation curve, and is restricted to a certain type of equation 
of state, herein we provide a general formalism that allows 
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us to deduce the density and pressure profiles without any 
prior assumptions about their shape or the equation of state. 

Analytic galaxy halo models that predict a signifi¬ 
cant amount of pressure or tension in the dark matter 
fluid include “string fluid” (Soleng 1995), or some varia¬ 
tions of scalar field dark matter (SFDM). See for instance 
Schunck (1999), Matos, Guzman & Urena-Lopez (2000), 
Matos, Guzman & Nunez (2000), Peebles (2000), and Ar- 
bey, Lesgourgues & Salati (2003). Our method provides a 
means of observing, or at least constraining, the pressure 
distribution in a galactic halo. Therefore it is in principle 
able to give evidence for or against specific proposed dark 
matter candidates. 

The key point is that in general relativity, density and 
pressure both contribute to generating the gravitational field 
separately. Furthermore, the perception of this gravitational 
field depends on the velocity of probe particles. These ef¬ 
fects become especially important when one compares ro¬ 
tation curve and gravitational lensing measurements, where 
the probe particles are fundamentally different: interstellar 
gas or stars at subluminal velocities for rotation curves, and 
photons which travel at the speed of light for lensing mea¬ 
surements. Our formalism accounts for these crucial differ¬ 
ences between the probe particles, and relates observations 
of both kinds to the the density and pressure profile of the 
host galaxy. Although we (mainly) consider static spheri¬ 
cally symmetric galaxies in a first post-Newtonian approx¬ 
imation, the basic concept is fundamental and can be ex- 
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tended to more general systems with less symmetry. A suit¬ 
able framework for considering most exotic weak gravity sce¬ 
narios is provided by the effective refractive index tensor, as 
introduced by Boonserm, Cattoen, Faber, Visser & Wein- 
furtner (2005). 

The present approach might also help to shed some light 
on prevailing problems that arise when combining rotation 
curve and lensing observations. For example, an unresolved 
issue exists when measuring the Hubble constant from the 
time delay between gravitationally lensed images: Using the 
standard models for matter distribution in the lens galaxy, 
the resulting Hubble constant is either too low compared to 
its value from other observations, or the dark matter halo 
must be excluded from the galaxy model to obtain the com¬ 
monly accepted value of Hq (Kochanek & Schechter 2004). A 
possible explanation of this trend might he in a disregarded 
pressure component of the dark matter halo. 

We organise this article in the following manner: First 
we introduce the minimal necessary framework of general 
relativity concepts, and point out the important condi¬ 
tions required to obtain the Newtonian gravity limit. Next, 
we elaborate on the post-Newtonian extension of the cur¬ 
rently employed rotation curve and gravitational lensing for¬ 
malisms. Consequently, we show how to combine rotation 
curve and lensing measurements to make inferences about 
the density and pressure profile of the observed galaxy. We 
then examine how noticeable the effects of non-negligible 
pressure could be in the measurements. Lastly, we dis¬ 
cuss how the formalism adapts to non-spherically symmetric 
galaxies and comment on the current observational situation 
and issues arising with the new formalism. 


2 GENERAL RELATIVITY FRAMEWORK 


In general relativity the motion of a probe particle is given as 
the geodesic of a curved space-time whose curvature is gen¬ 
erated by matter or more generally speaking, stress-energy. 
The static and approximately spherically symmetric grav¬ 
itational field of a galaxy is represented by the space-time 
metric (Misner, Thorne & Wheeler 1973, §23.2) 

ds 2 = — e 2Mr) At 2 + -- d? + r 2 dH 2 , (1) 

1 — 2 m(r)/r 

which is completely determined by the two metric functions 
<t>(r) and m(r). These coordinates (t,r,6,ip) are called cur¬ 
vature coordinates and 

dfi 2 = d(9 2 + sin 2 # dy> 2 (2) 


represents the geometry of a unit sphere. Invoking the Ein¬ 
stein field equations with the most general static and spher¬ 
ically symmetric stress-energy tensor gives the relation be¬ 
tween the metric functions and the density and pressure 
profiles: 


8n p(r) 
8np r (r) 
8np t (r) 


2 m'(r) 


i(r) 


_ r $'( r ) ( 1 _ 


2 m(r) 


— [m'(r) r — m(r)~\ [l + r <f?'(r)] + 


( 3 ) 

( 4 ) 
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+ F(r) 2 + $"(r) , 


( 5 ) 


where p(r) is the energy-density profile, and p r (r) and pt.(r) 
denote the profiles of the principal pressures in the radial 
and transverse directions. Note that we use geometrical units 
(c = 1 , Gn = 1) unless otherwise mentioned. Hence, if the 
metric functions $(r) and m(r) are given by observations, 
one can infer the density and pressure profiles. For a perfect 
fluid, one would expect p = p r = pt■ From (3) the physi¬ 
cal interpretation of m(r) as the total mass-energy within a 
sphere of radius r becomes clear. 


2.1 The Newtonian limit 

Standard Newtonian physics is obtained in the limit of gen¬ 
eral relativity where (Misner et al. 1973, §17.4): 

(i) the gravitational field is weak (2 m/r <C 1, 2$ <C 1), 

(ii) the probe particle speeds involved are slow compared 
to the speed of light and 

(iii) the pressures and matter fluxes are small compared 
to the mass-energy density. 

While there is no doubt that in a galaxy condition (i) is 
satisfied everywhere apart from the central region (Schodel 
et al. 2002), condition (ii) only holds for rotation curves and 
not for gravitational lensing. Finally, condition (iii) is related 
to the open question about the fundamental nature of dark 
matter. Hence, the possibility of dark matter being a high 
pressure fluid, or some sort of unknown field with high field 
tensions, cannot be excluded a priori. 

It is a standard result that condition (i) is enough to 
deduce that the gravitational potential <f>(r) is generated by 
the tt-component of the Ricci tensor (Misner et al. 1973): 

V 2 $ » R tt , (6) 

which on invoking the Einstein equations for R tt becomes 

V 2 4> « 4n(p + p r + 2p t ) . (7) 

Therefore, the metric function <f?(r) can be interpreted as 
the Newtonian gravitational potential 4 >n if and only if the 
pressures of the galactic fluid are negligible, i.e. if condition 
(iii) holds: 

V 2| I > n = 471 p. (8) 

It is now quite obvious from (7) that the gravitational field 
is highly sensitive to the pressure if density and pressure are 
of the same order of magnitude. 


3 ROTATION CURVES 


For the regime of rotation curve measurements, both condi¬ 
tions (i) and (ii) apply. In this case, the geodesic equations 
of the metric (1) reduce to (Misner et al. 1973) 


d 2 r 

dP" 


« -V4>, 


0 ) 


where r denotes the position vector of a probe particle. 
Equation (9) is equivalent to the Newtonian formulation of 
gravity, except for the general relativistic potential $ which 
replaces the Newtonian potential 4 >n- 
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Measurements of rotation curves are carried out by ob¬ 
serving the Doppler shift in the emission lines of the light 
emitting probe particles. In a general relativistic context, 
the observed shift in wavelength is not exclusively due to 
Doppler effects of the moving probe particles, but also de¬ 
pends on the gravitational redshift which arises as the pho¬ 
tons climb out of the gravitational potential well. 

It has been shown for edge-on galaxies that the total 
wavelength shift 1 z±(r) of an emission line of a probe parti¬ 
cle at radius r is given by (Nucamendi, Salgado & Sudarsky 
2001; Lake 2004; Faber 2006) 


1 + z±(r) = 


±|6| yjr 4>'(r) 


sjl — r \ e*W 


( 10 ) 


where prime denotes the derivative with respect to r, ' = 
d/dr and b is the impact parameter. z+ is the wavelength 
shift of an approaching particle and 2 _ that of a receding 
particle. 

The impact parameter is equivalent to the apparent dis¬ 
tance between the galactic centre and the emitting particle, 
once one takes notice of light bending effects. However, in 
the weak gravity regime of galaxies, where flat space is a 
suitable approximation, one finds 


|6| = r + (D[<S>\ 


( 11 ) 


for particles whose position vector r (with respect to the 
galactic centre) is perpendicular to the observer’s line of 
sight (Lake 2004). Thus, with the additional weak field as¬ 
sumption r<E>'(r) <C 1, equation (10) can be written as 


1 + z ± (r) = 1 =f \fr $'(r) + G[<I>, r 4>'] , (12) 

or equivalently, 

z± = r§\r) + 0[$> 2 ,(r$>'f /2 ,$\fr r ¥\ . (13) 


Comparing this expression to the Doppler shift in Newtonian 
gravity, 


= zf = r$' N (r), 


(14) 


we conclude that for small particle speeds v <C c, i.e. con¬ 
dition (ii), the observation of z± in edge-on galaxies is in 
first order equivalent to the Doppler redshift in Newtonian 
gravity, when $n is substituted by $. This also justifies the 
previous assumption r<f >, (r) -C 1. 

For galaxies of arbitrary orientation it is more tedious 
to obtain this result, but in a similar fashion, it can also be 
shown that the Doppler shift in wavelength is the dominant 
contribution to the observed total redshift (Faber 2006). 

Therefore, the usual techniques for obtaining the poten¬ 
tial Trc from rotation curve measurements can be employed 
if one keeps in mind that the motion of the observed parti¬ 
cles is not governed by the Newtonian gravitational potential 
<I>n, but by its general relativistic generalisation <E>: 


4>rc = 4> ^ $n • (15) 

If one assumes condition (iii), the density p is related to 4>rc 
by (8). In the general case, however, the interpretation of 


1 The wavelength shift that arises from the systemic velocity of 
the galaxy is not considered here, and Faber (2006) has shown 
that this does not change the result presented in this context. 


the mass which is inferred by rotation curve measurements, 
m-Rc(r), can be obtained from (7): 

w-rc (r) = r 2 $r C « 4n J (p + p r +2p t )r 2 dr. (16) 

Therefore, in the general case, we call mRc(r) ^ m(r) the 
pseudo-mass determined by rotation curve measurements. 


4 GRAVITATIONAL LENSING 

A fundamentally different approach of measuring the grav¬ 
itational field of a galaxy is gravitational lensing. Here, the 
observable photons are not only conveying the information 
about the gravitational field to us, they also act as probe par¬ 
ticles themselves. Hence, condition (ii) is naturally not sat¬ 
isfied for gravitational lensing observations. Consequently, 
the equations of motion for photons do not simplify to (9), 
as is the case for rotation curves. Instead, the geodesic equa¬ 
tions for photons have to be solved exactly to understand 
the influence of the gravitational field, as it is described by 
both metric functions, <f>(r) and m(r). Fortunately, for cer¬ 
tain spacetimes, such as e.g. (1), it is possible to characterise 
the entire trajectory of light rays with a single effective re¬ 
fractive index n(r). 


4.1 Fermat’s principle and the effective refractive 
index 

Fermat’s principle of shortest optical paths also applies to 
the geodesic trajectories of 4-dimensional curved spacetime 
(Kline & Kay 1965; Misner et al. 1973). This description of 
light rays in a gravitational field is equivalent to classical 
optics in a transparent medium with a continuous refractive 
index n, where Fermat’s principle is formulated as the van¬ 
ishing of the first variation of the optical length between two 
points, q\ and q 2 , on the trajectory: 

/■'32 

S n(f) [dr 2 + r 2 dff 2 ] =0. (17) 

Jqi 

By transforming the curvature coordinates of the spacetime 
(1) to so called isotropic coordinates (Perlick 2004), 

ds 2 = e 24,M {— dt 2 + n(r) 2 [dr 2 + r 2 dfl 2 ] } , (18) 

we introduce the scalar effective refractive index n(r) of a 
static spherically symmetric gravitational field. By direct 
comparison of (1) and (18), we find a differential equation 
that relates the r-coordinate of the isotropic coordinates to 
the r-coordinate of the curvature coordinates, 

dr f 

dr r\J\ — 2 m(r)/r 

and the refractive index, 

n[r) = — e . 
f 

Since condition (i) is satisfied for the region we are interested 
in, we can Taylor expand and formally integrate (19) under 
appropriate boundary conditions and find 



(19) 

( 20 ) 


( 21 ) 
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which, inserted into (20), gives 


n(r) 


exp 


<h[r(f)] 


m[r(f)l dr 

—--— — dr 

r(r) 2 dr 



( 22 ) 


where r(f) is given by the inverse of (21). Since f = 
r + 0[2m/r\, the radii in both sets of coordinates are in¬ 
terchangeable to the desired order and hence, we can also 
give the refractive index as a function of the curvature co¬ 
ordinate r directly: 


n(r) = 1 — <f>(r) 


f m{r) 

/ -Ar-dr + O 

( 2m \ 2 2m _. , 2 

- ,-$,4> 2 

] r 2 

V r ) r 


.(23) 


This effective refractive index entirely determines the tra¬ 
jectory of a light ray, i.e. the probe particles of gravitational 
lensing. Hence, it is the only possible observable of gravi¬ 
tational lensing. We note that the refractive index contains 
two distinct ingredients, the potential part, <f?(r), and the 
integral over the mass-function, J 2m(r)/r 2 dr. 

At this point, we conclude that since gravitational lens¬ 
ing observations yield n(r) and rotation curve measurements 
yield <f>(r), combined observations of n(r) and <f>(r) allow 
the separate deduction of <f>(r) and m(r), and therefore de¬ 
scribe the gravitational field of a galaxy in a general relativis¬ 
tic sense, without any prior assumptions. The fundamental 
principle is that the perception of the gravitational field by 
probe particles depends on the speed of the probe parti¬ 
cles, which manifests itself in the difference of observables 
n(r) ^ $(r). 

For convenience and comparability, we define the lens¬ 
ing potential as 


2 < f>iens(r) = 4>(r) + J m ^ 2 ^ dr , 


(24) 


so that 


n(r) = 1-2 <f>lens(r) + O [<f>? e ns • 


(25) 


4.2 Gravitational lensing formalism 


The standard formalism of gravitational lensing in weak 
gravitational fields is based on the superposition of the de¬ 
flection angles of many infinitesimal point masses (Schnei¬ 
der, Ehlers & Falco 1992, §4.3). 

In general relativity, a point mass M is described by the 
Schwarzschild exterior metric, 


ds 2 



dt 2 + 


dr 2 


1 - 2 M/r 


+ r 2 dff 2 . 


(26) 


which inserted into (23) gives the effective refractive index 


whereas from (27) it is clear that in the general case, the 
refractive index is only partially specified by the potential 
term. That the mass term of the refractive index of a point 
mass is identical to the potential term is a special case of 
the Schwarzschild metric. Keeping this in mind, we proceed 
to outline the current formalism. 

For a point mass, the angular displacement of the grav¬ 
itationally lensed image in the lens plane can be calculated 
from the refractive index (28) (Schneider et al. 1992): 

« = 4M ’ (3°) 

where £ is the vector that connects the lensed image and the 
centre of the lens in the 2-dinrensional lens plane. Since the 
extent of a lensing galaxy’s mass distribution is small com¬ 
pared to the distance between the light emitting background 
object and the lens, as well as compared to the distance be¬ 
tween the lensing galaxy and the observer, one assumes the 
deflecting mass distribution to be geometrically thin. There¬ 
fore, the volume density p of the lensing galaxy can be pro¬ 
jected onto the so-called lens plane, resulting in the surface 
density E(£) which describes the mass distribution within 
the lens plane (Schneider 1985). 

The total deflection angle of a lensing mass with finite 
extent is then said to be given by the superposition of all 
small angles (30) due to the infinitesimal masses in the lens 
plane (Schneider et al. 1992): 

&(o=4 yW) ^~^ |2 dv. (3i) 

This equation is the foundation of the gravitational lensing 
formalism as introduced by Schneider (1985), and Blandford 
& Narayan (1986). 

Since this formalism is based on the assumption that 
the total angle of deflection is caused by the superposition 
of point masses - without the notion of pressure at all it 
automatically assumes that the underlying lensing potential 
is Newtonian, i.e. <f?i ens = 4>n. Hence, the lensing potential 
and the naively inferred density- or mass-distribution are 
related by (8): 


^ ^ > lens(^’) — 47tpiens(^) 

(32) 

which implies 



(33) 


However, as we argued previously, the lensing potential ‘Fiens 
is the fundamental observable, and not the density p which 
was used to construct the formalism. Therefore, for the gen¬ 
eral case that does not assume (iii), i.e. Tiens ^ 4?n, we note 
that 

piens(r) ± p{r) and miens(r) ^ m(r). (34) 



In the Newtonian limit, this is generally identified with the 
Newtonian potential, 

n(r) = 1 — 2 <f> N (r), (29) 


Instead, the deduced mass distribution mi ens (r) has to be 
considered as a pseudo-mass similar to that of rotation curve 
measurements. Its physical interpretation can be deduced 
from the definition of the lensing potential (24): 


miens (r) = i mnc(r) + i m(r) 


471 


P+ \ (Pr + tyt) 


(35) 

(36) 



5 BRINGING ROTATION CURVES AND 
GRAVITATIONAL LENSING TOGETHER 


We showed in the previous sections that the potentials ob¬ 
tained from rotation curve and lensing observations, 3 >rc 
and 'hiens, do not agree in the general case, 


< I > Rc(r) 

= $(0> 

(37) 

4>lens(r) 

1 . . if m(r) 

= 2 $ «+2 / r 2 

(38) 


but only in the Newtonian limit, where condition (iii) holds, 
in which case 3 >rc = 'f’lens = $n- Since this is the standard 
assumption for interpreting rotation curve and lensing data, 
the results of these observations are often reported as mass 
distributions instead of potentials. Under the Newtonian as¬ 
sumption, the mass and the potential are related by a field 
equation of the form (8). In the general case, this leads to 
the definition of the distinct pseudo-masses, 


mRc(r) = r 2 <f>'(r), 

(39) 

1 2 / 1 

m iens (r) = - r $ (r) + - m(r), 

(40) 


which describe the observations equivalently to the poten¬ 
tials (37) and (38). Equations (39) and (40) can easily be 
inverted to give the metric functions 4>'(r) and m(r), 


s mRc(r) 

$ (r) = J , 

(41) 

m(r) = 2miens(r) - m R c(r), 

(42) 


which inserted into the field equations of general relativity 
(3)-(5) yield the density and pressure profiles: 


4 nr 2 p(r) = 2mL(r)-ml lc (r), 
4nr 2 p r (r) = 2 ™nc(r) - m lens (r) + p 


(^y 


47t rpt (r) = r 


nine (r) - mi ens (r) 


/ 

+ o 

/2 m\ 2 


\ r ) 


= [4tt r 2 p r {r)]' + 0 (/^) 
As consistency checks, we note that: 


(43) 

(44) 


(45) 


• The Einstein equations in curvature coordinates, (43)- 
(45), agree to the given order of 2 m/r with the Einstein 
equations of the metric in isotropic coordinates, and there¬ 
fore the approximation r w r is valid; 

• From the equations (43)-(45) follows that 

47t r 2 \p(r) +Pr(r) + 2 pt(r)] « m' RC {r) , (46) 


and 


4tt r 


p(.r) + - (p r (r) + 2p t (r)) 


m[ ens (r ). 


(47) 


and thus these results are consistent with the weak field 
approximation of the field equations (7), and the interpre¬ 
tations of the pseudo-masses (16) and (36); 

• For m.R C (r) = m] ens (r) = m'(r ) we find the desired 
result of the Newtonian limit: 


47t r 2 p(r) 
4nr 2 p r (r) 


m'(r), 



(48) 

(49) 
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Figure 1. The x-f&ctor (53) as a function of w in the commonly 
discussed range w G [—1,1]. Naturally x(0) = 1, which corre¬ 
sponds to the Newtonian case, see (48) - (50). At the ends of the 
plotted range we find x(—1) = 1/4 and x(l) = 5/8. The zero¬ 
crossing and the first-order pole correspond to w = —2/3 and 
w = —1/3 respectively, as is obvious from (53). 


47tr 2 p*(r) 



(50) 


We conclude that the currently existing formalisms for 
analysing data from rotation curve and gravitational lensing 
observations can be used to separately obtain the pseudo¬ 
masses mRc(r) and miens (r), which by (43)-(45) yield the 
density and pressure profiles in a first post-Newtonian ap¬ 
proximation. Furthermore, from the combination 

47tr 2 (p r + 2p t ) « 2 (m' RC - m[ ens ), (51) 


one can immediately infer that the observed system is New¬ 
tonian in the sense of condition (iii) if and only if rn' nc (r) ~ 
m[ ens (r). Furthermore, defining the dimensionless quantity 


w(r) = 


p r (r) + 2p t (r) 2 mR C (r) - m( ens (r) 


3 p(r) 


3 2 Sens ( r ) - «4c( r ) 


(52) 


gives a convenient parameter that determines a “measure” 
of the equation of state. 


6 PARAMETERIZING THE SIZE OF THE 
EFFECT 


To get an idea how noticeable the existence of a pressure con¬ 
tribution is likely to be in the measured data, we introduce 
the X"f ac tor which we define as the ratio of the derivatives 
of miens and mRc: 


xHS 


SensS 

SicM 


2 + 3ui(r) 
2 + 6w(r) 


(53) 


This can easily be obtained from rearranging (52). 

Using fig. 1, we can now see how the ratio of the slopes 
of m lens and m rc relates to w. One should especially note 
that x i s n °t very different from unity in the vicinity of 
w = 0, making it difficult to detect small pressures. However, 
for w € [1/2, 1], which is a range that might plausibly be 
identified with real-life data (or at the very least with some 
real-world theoretical prejudices), x takes values between 
62.5% and 70%. For most negative w, specifically for w in the 
interval [—1,0], x is rather distinct from unity, and therefore 
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should be easily detectable if the quality of the observational 
data is high. 

Also note that in the special case of a position- 
independent w — so that the “equation of state of the galac¬ 
tic fluid” is constant throughout the observed region — the 
ratio between mi en s and to-rc will be the same as between 
their derivatives: 

miens(r) = x(w) WRc(r) • (54) 

This relation is likely to be more useful since it does not 
depend on numerically obtained derivatives when one wishes 
to compare mass profiles. However, it comes at the price of 
the additional assumption that w = constant. 


7 NON-SPHERICAL GALAXIES 

Although we chose to present this new formalism using the 
example of a simple spherically symmetric galaxy, it is easy 
to show that the fundamental concept behind the formalism 
is also valid for configurations with less symmetry. 


7.1 Rotation curves 


The formalism of Newtonian mechanics is usually adopted 
when data from dynamical observations is examined to de¬ 
termine the shape of the matter distribution of a galaxy. 
That is, the fundamental equation employed is (9): 


d 2 r 

d IF 


-V${x,y,z). 


(55) 


Note that no particular symmetry is now assumed for the 
gravitational potential $(x,y,z). The only necessary as¬ 
sumptions are the weakness of the gravitational field [con¬ 
dition (i)] and the slowness of the probe particles [condition 

(ii)l- 

In view of the formalism presented herein, one now only 
needs to realise that the gravitational potential inferred from 
dynamical observations is generally not the Newtonian po¬ 
tential 


$rc (x,y,z) = ${x,y,z) ^ $ N {x,y,z ), 
and that the valid field equation is [cf (7)]: 


V 2< h(x, y , z) « 47t 


p(x,y,z) +2^._ Pi{x,y,z) 


(56) 


(57) 


The first Einstein field equation 2 is then 

4 np{x,y,z) = -V 2 <I >{x,y,z) - V 2 h(x,y,z) 

+0[$ 2 ,h$,h 2 ] . (61) 

Inverting this equation yields the effective refractive index 

n(x, y, z) = 1 - <b{x,y,z) - 4tt (V 2 ) _1 p(a:,y, 2 ) 

+0 [<f? 2 , /i 2 ] , (62) 

where the constants of integration have been chosen to 
agree with the special case of spherical symmetry (23), and 
(V 2 ) -1 is the inverse Laplacian operator. The general non- 
spherical lensing potential analogous to (24) can be defined 
as 


so that the non-spherical refractive index is of the same form 
as (25): 

n{x,y,z) = 1 — 2 <l?i ens (a;, y,z) + 0 [$Ls] • (64) 

corresponding field equation for the lensing potential is 

1 V7 2 ,T./~ „. -V I ~ ... g\ (65) 


The 


V 2 $! 
= 471 


'lens (x,y,z) = ^V 2 $(x,y,z) + 2np(x,y, 

r 1 v^ 3 n 

p(x,y,z) + - > Pi(x,y,z) 

Z zJ i=l 


( 66 ) 


7.3 Non-spherical formalism 

Combining the non-spherical field equations (57) and (66) 
yields the density and pressure distributions in absence of 
any particular spatial symmetry: 

4np « 2 V 2, I>ie n s(a:, y, z) — V 2< f>Rc(a:, y, z) (67) 

47t pi « 2 [V 2 <J>Rc(a:,y, 2 ) - V 2 4>i en s(a:, y,«)] (68) 

Thus the two observable potentials 4 >rc (x,y,z) and 
4 > iens(*, y, z) determine the density and pressure distribu¬ 
tions of a non-spherical galaxy to the lowest order in 
the weak gravitational field represented by the functions 
&(x,y,z) and h(x,y,z). This is the straightforward exten¬ 
sion of the spherically symmetric formalism presented in this 
paper. 


8 OBSERVATIONAL SITUATION 


7.2 Gravitational lensing 

The notion of a simple effective refractive index can easily be 
extended to the class of conformally static spacetimes which 
are also conformally Euclidean (Perlick 2004). The metric of 
this class has no particular spatial symmetry and takes the 
form: 

ds 2 = e 2 ^ x ’ v ’ z l {—dt 2 + n(x,y,z) 2 [da; 2 + Ay 2 + dz 2 ] } .(58) 

To account for the weakness of the gravitational field [con¬ 
dition (i)[, we assume 

$(*,y,2)< i (59) 

and 

n(x, y, z) = 1 + h(x, y, z) with h(x, y, z) <g 1. (60) 


The post-Newtonian formalism we have outlined requires 
the simultaneous measurement of (pseudo-)density profiles 
from rotation curve and gravitational lensing observations. 

While in principle these profiles do not have to be of the 
same galaxy, they must be comparable in the sense that they 
accurately describe “similar” galaxies. For example, weak 
lensing measurements can be used to statistically infer the 
(pseudo)-density profile of an “average” galaxy (Brainerd 
2004). At the same time, analysing the dynamics of satel¬ 
lite galaxies gives the rotation curve and thus, the corre¬ 
sponding pseudo-density profile, of another “average” galaxy 
(Brainerd 2004). Whether these two “average” galaxies are 
comparable or not depends on many factors, such as e.g. 

2 The first Einstein field equation is that which is associated with 
the /./-component of the Einstein tensor. 
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the distribution of galaxy morphologies in both samples, 
the statistical noise, the employed models for the (pseudo)- 
density distribution, etc. These statistical issues render the 
fast-growing collection of weak lensing data problematic for 
our purposes. 

On the other hand, combined simultaneous measure¬ 
ments of rotation curves and lensing of individual galaxies 
are extremely well suited for our formalism. However, while 
there is a large number of individual rotation curves avail¬ 
able (> 100,000; Sofue & Rubin 2001), the number of in¬ 
dividual “strong” lensing systems with multiple images is 
rather limited 3 (~ 70; Kochanek et al. 2005). Combined ob¬ 
servations are further aggravated by the differing distance 
scales: Most high quality rotation curves are naturally avail¬ 
able for galaxies with a low to intermediate redshift of up to 
z ~ 0.4 (Sofue & Rubin 2001), while gravitational lenses are 
easier to detect at intermediate to high redshifts (z > 0.4; 
Kochanek et al. 2005), since the image separation scales in¬ 
creasingly with the redshift of the lensing galaxy (Schnei¬ 
der 1985; Kochanek & Schechter 2004). Therefore, even for 
nearby galaxies with existing combined measurements of 
kinematics and lensing (e.g. 2237+0305 at z ~ 0.039 and 
ESO 325-G004 at 2 « 0.035), the lensing data is restricted 
to the core region, while the rotation curve is only described 
by few data points in the outer region of the lens galaxy 
(Barnes et al. 1999; Trott & Webster 2002; Smith et al. 
2005). Consequently, the inferred pseudo-mass profiles are 
available for the same galaxy, but unfortunately at different 
radii and therefore not comparable. 

Although the observational situation makes it currently 
difficult to employ the formalism presented, the situation is 
likely to improve in the future when observations with a 
higher resolution will be carried out - preferably with an 
emphasis on obtaining high-resolution rotation curves for 
lensing galaxies that exhibit lensed images at different radii. 


9 CONCLUSIONS 

We have argued that the standard formalism of rotation 
curve measurements and gravitational lensing make an a 
priori Newtonian assumption that is based on the CDM 
paradigm. We introduce a post-Newtonian formalism that 
does not rely on such an assumption, and furthermore allows 
one to deduce the density- and pressure-profiles in a gen¬ 
eral relativistic framework. In this framework, rotation curve 
measurements provide a pseudo-mass profile m.Rc(r) and 
gravitational lensing observations yield a different pseudo¬ 
mass profile miens (r). Combining both pseudo-masses allows 
one to draw conclusions about the density- and pressure pro- 


files 4 in the lensing galaxy, 



p(r) 

1 

47tr 2 

[2 m[ en s {r ) — m,R C ( r )] , 

(69) 

p r (r) + 2 p t (r) 

2 c 2 
47tr 2 

["4c (r) — m[ ens (r)] . 

(70) 


In the case of absent or negligible pressure, this could be used 
to observationally confirm the CDM paradigm of a pres¬ 


3 For an up to date list see 

http://www.cfa.harvard.edu/glensdata/. 

4 These formulae are given in SI units, hence the factor of c 2 . 


sureless galactic fluid. Conversely, if significant pressure is 
detected, a decomposition of the galaxy morphology would 
allow new insight into the equation of state of dark matter. 

For instance, detailed observation of the recently dis¬ 
covered closest known strong lensing galaxy ESO 325-G004 
(Smith et al. 2005) could provide satisfactory data to allow 
the decomposition of density and pressure of the galactic 
fluid, as outlined in this article. The system consists of an 
isolated lensing galaxy at redshift 2 « 0.035 with an effec¬ 
tive radius of R e g = 12'/5 and arc-shaped images of the 
background object at R ss 3”, and possible arc candidates 
at R ~ 9”. Smith et al. (2005) intend to collect more de¬ 
tailed data that hopefully will include extended stellar dy¬ 
namics and hence, allow for a direct comparison of the rota¬ 
tion curve and lensing data, if the arc candidates at R « 9” 
turn out to contribute to the measurements. 

Since the formalism presented is based on a first-order 
weak field approximation, we suggest that to confirm the 
findings, one should re-insert the obtained density and pres¬ 
sure profiles into the metric (1). The actual observed quanti¬ 
ties can then be extracted numerically for comparison from 
the exact field equations (3)-(5) and the geodesic equations. 

Finally, even though data might not yet be available to 
constrain the dark matter equation of state noticeably, one 
should note that the possibility of non-negligible pressure 
in the galactic fluid introduces a new free parameter into 
the analysis of combined rotation curve and lensing obser¬ 
vations. 
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